Electrons on a sphere in disorder potential 
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We investigate, both analytically and numerically, the behavior of the electron gas on a sphere 
in the presence of point-like impurities. We find a criterion when the disorder can be regarded as 
small one and the main effect is the broadening of rotational multiplets. In the latter regime the 
statistics of one impurity-induced band is studied numerically. The energy level spacing distribution 
function follows the law P(s) ~ sexp(— as h ) with 1 < b < 2. The number variance shows various 
possibilities, strongly dependent on the chosen model of disorder. 
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I. INTRODUCTION 

A considerable amount of the theoretical interest has 
been devoted last years to the investigation of electronic 
properties of the nano-size objects possessing regular 
cylindrical and spherical shapes. While most of these 
studies were consecrated to the various properties of car- 
bon nanotubes, the technological advances in production 
of spherical objects also require further theoretical under- 
standing of peculiar quantum effects of electronic motion 
in such materials. 

A well-known fullerene molecule, Cgo, could be one 
of the examples, along with its further modification, 
so-called "onion" graphitic structures. Q The spherical 
nano-size objects are also found in the studies of the non- 
linear optical response in composite materials Q and 
of simple metal clusters. || A rapidly evolving field of 
photonic-band-gap materials [0] provides yet another ex- 
ample of spheres of about 300 nm in diameter. In the 
latter case, one may find silica balls with semiconducting 
coating (coated opals) ||, or "inverse" opal structures, 
where the initial SiC>2 template is chemically removed 
and carbon spherical shells form three-dimensional fee 
structure. 

In all these cases one can assume that the motion of 
itinerant electrons is confined within a spherical layer of 
a width small comparing to a radius. To a first approxi- 
mation, one ignores the interaction effects and considers 
a situation of the electron gas. Theoretical efforts in this 
direction comprise an analysis of the behavior of such gas 
on a sphere in the uniform magnetic field and an evalua- 
tion of the electronic correlations without a field. It was 
shown particularly that if the field is small so that the 
magnetic length is larger than the radius, then one can 
expect the jumps in the magnetization and susceptibility 
of the sphere. (?],|^] An extension of this analysis for the el- 
lipsoid of revolution has been undertaken recently. pJTo| 
In stronger fields, when the magnetic length is smaller 
than radius, one finds a series of interesting effects de- 
scribed elsewhere. Jll|,|| 

The electronic correlation functions for the topology 



of the sphere exhibit, particularly, a peculiar coherence 
effect when the coherence length of electronic motion ex- 
ceeds the radius of the sphere. It turns out, that 
the amplitude of these functions is enhanced for the an- 
tipodal points on the sphere, where all partial waves of 
quantum motion come in phase. 

The electron gas approximation, employed in these 
works, should be violated in experimental realizations of 
nanospheres. One of the reasons for this violation is the 
effects of electron-electron interaction, which were stud- 
ied in ||. 

Another source of possible inadequateness of the elec- 
tron gas model is the explicit absence of rotational sym- 
metry, inferred by the presence of impurities. Indeed, the 
experimental realization of the spherical layer, wherein 
the electrons are confined, might be far from the ideal 
shape. The inhomogeneities of various kind (impurities, 
fluctuations of layer's width) should eventually destroy 
the effects, found theoretically for the idealized model. 
As a first crucial step here, one observes that the im- 
purities necessarily raise the multiple degeneracy of the 
energy levels, found initially in a quantum rotator model 
of a free gas. 

In this paper we consider the effects of potential impu- 
rity scattering in the electron gas on a sphere. We em- 
ploy the model of point-like impurities, which is shown 
to be valid if the angular momentum of an electron does 
not exceed inverse angular range of the potential. Our 
results suggest that there is a definite range of param- 
eters, where the electron gas approximation should be 
applicable. Particularly, this approximation is validated, 
when the radius of the sphere is smaller than the ra- 
tio \pvfvi mv , with the areal densities of electrons and 
impurities, v and Vimp, respectively. We show that in 
this case the degeneracy of the initially degenerate mul- 
tiplet, close to Fermi energy, is lifted only partially. The 
unsplit states are superpositions of the initial spherical 
harmonics, and hence describe the electrons freely prop- 
agating along the sphere. As a result, one expects, that 
for sufficiently small disorder it is possible to discuss the 
coherence effects, induced by the spherical topology. 

For the regime of weak disorder, we numerically in- 
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vestigated the splitting of on multiplet, caused by ran- 
dom impurities. We found that the statistics of the 
energy levels in this case is not identical, albeit close, 
to the predictions of the random matrix theory for the 
Gaussian orthogonal ensemble. Particularly we show 
that the energy level spacing distribution follows the law 
P(s) ~ s exp — s b , with b ~ 1.76 in the limit of white- 
noise potential distribution. We discuss a possible reason 
for this deviation, related to the geometry of the problem. 

A rest of the paper is organized as follows. We formu- 
late the problem and treat the random potential pertur- 
bativcly in Sec. II. We identify here the region of param- 
eters, where the perturbation theory fails. In Sec. Ill we 
numerically investigate this region. The discussion and 
conclusions are presented in Sec. IV. 



II. PERTURB ATIVE TREATMENT OF THE 
IMPURITY POTENTIAL 



A. Green's function 



In the ideal situation the motion of an electron on the 
sphere is described by the quantum rotator model. The 
wave functions W are the spherical harmonics Yi m and 
the spectrum is simple : 



We normalize VP ( 



), Ei = {2m e rl)-H{l + 1). (1) 
on the surface of the sphere 



The main effect of the random impurity potential is 
the lifting of (21 + 1)— fold degeneracy of each I— th en- 
ergy level. Qualtitatively, one expects that at sufficiently 
small disorder these split levels form a subband nearly 
the initial position Ei, and the width of this subband 
does not exceed the separation between the adjacent lev- 
els Ei and Ei±\. This regime corresponds to the mean 
free path of the electron, l m fp, larger than r$. The co- 
herence effects of the electronic motion should be still 
present in this case. plJl^] 

At larger disorder, the width of induced subband is 
more than \Ei — Ei±\\. The different subbands over- 
lap now and, as a result, one finds a constant density 
of states, typical for a planar two-dimensional electron 
gas. Having the relation l m f p < r$ in this case, one sees 
that the spherical surface is effectively decomposed onto 
the "patches" of size Z m / P , in which the electronic motion 
is essentially planar. 

To find out the corresponding criterion for these two 
regimes, we use the Green function formalism. Treat- 
ing disorder as perturbation, we write the self-consistent 
equation for the Green function. The spread of the split 
levels around the position of the initial multiplet is de- 
scribed by an imaginary part of the corresponding self- 
energy. This approach is very similar to one employed 
in the analysis of splitting of the lowest Landau level 



in two-dimensional electron gas (2DEG), [jl5 14 and a 
reader can easily find parallels with our situation below. 
The Green function is given by 



*r m (0',<«*z m (0,o) 



U3 + [I - E lm 

For the free motion this expression simplifies 

2 ^_ lY ^ (21 + l)P z (cosft) 
to + (j, — Ei 



G°H = (4 7 rr 2 )- 1 X: 



(2) 



(3) 



with the distance between two points on the sphere 
cos ft — cos 9 cos 8' + sin 9 sin 8' cos cf>. 
The basic equation is 

G(0 12 ) = G°(0 12 ) + y"rfr 3 rfr 4 G (r!i 3 )S(r!34)G(r!42) 

(4) 

where we write f2i2 for the angular distance between Qi 
and ^2, etc. The integration over the surface of the 
sphere reads as j dr = r$ J dfl. 

Being averaged over disorder, the functions entering 
the above equation depend only on the distances between 
the corresponding points. Representing each function 
through its generalized Fourier coefficients, 



F(Q) 



2^ + 1 
An 



FiPi(Q), 



we have 



Gr 1 = (G?)- 1 - 



(•5) 



(6) 



with = (u> — Ei)^ 1 . We consider different contribu- 
tions to Tii below. But before doing that, we list some 
formulas applicable to the spherical geometry in the next 
subsection. 



B. some useful formulas 

First we note that eq. allows other representations, 
namely 



F(fii 2 ) = r - 2 ^F ; A ; y ;o (f! 1 2), 
/ 

= r - 2 ^F ; y ; ; l (r! 1 )y ;m (r! 2 ), 



In i 



with the quantity 

Ai = y/(2l + l)/(4n). 



(7) 



(8) 



The Green function (^) can be represented || through 
the Legendre function 



G°M = - 



2 cos7rA 



P_i/ 2 +a(-cosQ), 



(9) 
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with A = y2mero(/"t + w) + 1/4. In subsequent consider- 
ation we are mostly interested in energies close to Fermi 
level, u> <C /i. In this sense, we define the angular Fermi 
momentum, L, as -El — A*, more precisely 



A ~ L + 1/2+ (w - Ei)/AE 
AE = E L+1 -E L = L/{m e rl). 

For small angular distances ft -C A^ 1 we have 

ft 



me 
2tt 



2 In A sin ■ 



7 — 7r tan 7r A 



(10) 
(11) 



(12) 



with the Eulcr constant 7 = 0.577 . . .. This expression is 
similar to the Green function for the 2DEG in the mag- 
netic field B. Indeed, in the latter case we have for r -C I* 

M 



to 

7r tan7r — 



(13) 



where the magnetic length = (eS) -1 ' 2 and the cy- 
clotron frequency lu c = eB/m e . 

For larger distances, A sin ft > 1 we have approxi- 
mately 



G 



cos[A(tt - ft) - tt/4] 



V27rA sin ft 



cos 7rA 



(14) 



As was shown in ||, the existence of two oscillating ex- 
ponents in (JOJ) , exp±iAfti2, indicates the interference 
of the partial waves on the sphere, one wave going along 
the shortest way between the the points fti and ft2, and 
another going along the longest way, turning around the 
sphere. 

We discuss now a notion of point-like impurity on 
the sphere. Letting first the range of the potential to 
be zero, we write uo8(r — r^) = UQrQ 2 S(fl — ftj) = 
uorQ 2 S(cos8 — cos6i)5(4> — 4>i) . Next we introduce the 
function describing the formfactor of impurity in the form 
<5 a (ft) = (27ra)-V cose - 1 )/ a , the limit a < 1 is assumed. 
The physical range of the potential is rQ^/a. To exponen- 
tial accuracy we have J dflS a (fl) — 1. Then, using the 
asymptotic (Macdonald) formula for the Legendre poly- 
nomials at small 9, we calculate the auxiliary integral 



A, = Al 1 / dftr /0 (ft)(5 a (ft) ~ exp 



1(1 + 1) 



(15) 



This equation shows that as long as considered moments 
/ < 1/i/a, one can approximate the potential by delta 
function. 

The impurity potential centered at the north pole reads 

as 

u(ft) = u t r^ 2 S a (n) = u r^ 2 ^ A z A i r i0 (ft) (16) 



and the potential centered at the point ft^ is 



t (ft - ft,) = mr^Y, A i r i ; i (ft l )r /m (ft). (17) 



The averaging over the disorder (. . .) is written as 



Zi(ftl l )u(ft2i) = Vim p r / (2ftj1t(ftli)lt(ft2i 



- 2 ^A 2 A ; r ;o (fti 2 ) 



r <5 2a (fti2) 



(18) 
(19) 

(20) 



Here Vi mp = Ni mp / (Airr 2 ,) is the impurity concentration 
and Ni m p is the total number of impurity centers. 

Consider next more complicated objects, which are 
to be used below. If we have two functions, /(ft) = 
Sim flmAiYi m (Q.) and g(ft) = Y^i m 9im,AiYi m (Q), then 



/(ft)g(ft) = ^ (f hmi 9l 2 m 2 V^Z\ m2 ) A ls Y l^M 

li.rrij 



V, 



— a 2 a 2 A- 2 r< l 3 r il ^ m 3 



limil27Ti2 h I2 I3 /1O/2O livciil^m 



(21) 



with the Clebsch-Gordan coefficients Of™ l2 . 0] 

Particularly, when /(ft) = S a (Q) and <?(ft) can be ap- 
proximated by a constant g at small ft < s/a, we have 
(/ ' g)lm = gAi5 m0 . 



C. perturbation theory 



We use the model of point-like impurities randomly 
distributed on the sphere. The amplitudes uq of the 
impurity potential are chosen to be equal. The overall 
potential is 

C/(ft) = u r- 2 A/^(ft) Y U^i) (22) 



im 



i=l 



As usual, the average of this potential gives a constant, 
which is incorporated into the chemical potential and is 
discarded below. First nontrivial diagrams in the pertur- 
bation series in uq are shown in Fig. [jj 



.«3 



a b c 

FIG. 1. First diagrams describing impurity scattering 
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In the self-consistent Born approximation, we find for 
the self-energy in the second order of U, Fig. |l|a 



so that e[°^ is independent of I (as long as A; ~ 1) 



,(a) 



i 2 G(n ~ Jo) 



(23) 



(24) 



The last line is obtained upon an assumption L^/a < 1, 
when the argument of the logarithm in ( |l2|) is small. Re- 
calling that L <-> kpro, the latter condition means that 
the range of the impurity potential does not exceed kp 1 
in real space. 

Letting lo = and introducing the quantity 

z = ttZi/AE 

we find the self-consistency equation in this order of U : 



: ~ Xu 2 ( — \n(L\/a) — cot z ) , 



(25) 



with 



X = N imp /(2L + l), (26) 

and the dimensionless scattering amplitude 

u = UQm e /2. (27) 

In the following we adopt the picture of the weak scat- 
tering, u<l. 

The real part of in (^5|) redefines the chemical po- 
tential and may be ignored. The imaginary part of 
describes the width of the impurity-induced band, aris- 
ing instead of degenerate multiplet L : 

z = Xu 2 {\ ln(LVa) - i) , Xu 2 > 1 (28) 

(29) 



= -iVXu, Xu 2 < 1 



As we discussed above, the second regime is of our main 
interest. It corresponds to the case when the cot z dom- 
inates over the logarithm in (|25|), which in turn means 
the possibility to ignore the transitions between the lev- 
els with different I. The scattering in this case is between 
the states with different m within one multiplet I. 
The next diagram, Fig. |l]b, is wrtten in the form 

£W(fi) ~ y lmp ulr^ 2 5 a {n)G 2 {n). (30) 

Or, in the previous notation 



7r£, (6) , /2 , ^ 

~ Xu 3 - ln(Ly/a) - cot ; 

AE \ 7T 



(31) 



In order to find the relative importance of this cor- 
rection, we put the former estimate, S a , into the energy 



argument of the Green function, A ~ I + 1/2 — Y,f/AE. 
Then we find 



Ef/E" ~ u, Xu 2 > 1, 

-X- 1 ' 2 , Xu 2 <l. 



(32) 
(33) 



Evidently, the correction E & can be neglected at large X 
and dominates at X < 1. We discuss this feature in more 
detail below. 

The next diagram is shown in the Fig. |l|c 



E c (tt) = v?utG{n) 



(34) 



The Fourier component E^ behaves differently in the two 
above cases z > 1 and z < 1. 

In the case z«l the Green function can be approxi- 
mated by the expression p). Near the level El we have 



G ~ (47rr 2 ) 



„ 1 (2X+l)P L (cosfi) 



lo + fi - E L - E L 
m e cos[(L + l/2)f7-vr/4] 
z y/2TrL sin £1 



(35) 

L sin ft > 1 (36) 



where we have put lo = in the last line. Using ( p5| ) and 
(|36|), one finds for T,f the following expressions, respec- 
tively, 



4 2 f 



'*n P »0'0 



2 imp "o^ 1 ; 



/ . v 7i.0i0 v i0i0> 



-AE ^^ln 



2ttzA l 
L 



In 



L-Z 



ZLz 3 



(37) 

(38) 
(39) 



Note that the expression (|37|) tells us that : i) I and L 
are even (odd) simultaneously, i.e. (L + 1) = 0mod2, and 
ii) E^ = at Z > 3L. The Eq. (|37| ) is not convenient 
and one can use an approximate expression (^). In the 
considered case (E9f) we estimate at I = L : 



(40) 



E^/E£ ~i- 1 lni< 1 



Therefore this correction is smaller than the previous one. 

The case Im z > 1 corresponds to the mean free path 
of the electron smaller than tq. In this case the Green 
function ( |l4| ) exhibits the presence of only one wave, with 
the shorter path. One has from (nj) at ImX > 1 : 



Q ~ e 

%/2^Af2 
As a result we have 



exp 



/(*) 



AE XW /< + l/2 
' J V 3A 
2x 



(2A) 3 
V6(l + ^)^ 



-.K 



1 



(41) 

(42) 
(43) 
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with the complete elliptic integral K (x) and A ~ L + 
1/2 + z/n. A comparative importance of this contribu- 
tion is 



(44) 



and may be significant in the considered case Xu 2 > 1. 

One can evaluate also a sixth-order correction, shown 
in the Fig. ^, which is the next non-trivial term. In case of 
our interest, Xu 2 < 1, the contribution of this diagram, 
£ d , is estimated as follows. 



AE- 



X A u 



3^6 



(45) 



so that the relative importance of this diagram is again 
small : 



££/£» ~ 1/L 



/ v \ ✓ \ 



(46) 



FIG. 2. A sixth-order diagram with maximally crossed im- 
purity lines 



One can further argue that the weak localization cor- 
rections do not contribute much at Xu 2 < 1. Indeed, in 
this case of our interest below, when the mean free path 
is larger than the radius, we are essentially in the ballistic 
regime, hence one does not find a typical enhancement 
in the cooperon series. As a result, the arguments eluci- 
dating the role of the spherical topology in the diffusive 
regime Q cannot be fully applied to our case. 

On the other hand, let us consider a diagram which 
describes four scattering events on the same impurity. 
This is the next diagram in a sequence, formed by Figs. 
|l|a, |l|b. It can be easily shown, that in the regime 129) 
this diagram is estimated as ££/A, being thus large at 
X < 1 . This latter sequence of diagrams contributes to 
the S-matrix for one impurity. The formal divergence 
in this sequence at X < 1 indicates a singularity in the 
S-matrix. In the following section we show, that this 
singularity corresponds to the incomplete raising of de- 
generacy in L— th multiplet. 

Note that the quantities L, X and (A_E) _1 scale as ro 
so that in the limit of large sphere the regime (P9|) docs 
not occur and we are left with eqs. (|28|), (p2|), (|44D. It is 
worth to represent the quantity X in the form 



X = 



'imp 



T , 



(47) 



then we see that the criterion X < 1 means ro < 
y/v/vimp- For the semiconducting situation, adopting 
a simple estimate v ~ Vimp ~ 10 10 cm -2 , we get for this 
regime ro ^ 100 nm, the latter value comparable to the 
size of opal balls. 



Concluding this section, we observe that at sufficiently 
small disorder, Xu 2 < 1, the impurity-induced bands, 
referring to different multiplets, do not overlap. This 
regime corresponds to the scattering mostly within one 
multiplet. At small X one meets a singularity in S-matrix 
of scattering for one impurity. We consider it to some de- 
tail in the next section. 



III. SCATTERING WITHIN A MULTIPLET 

In this section we restrict ourselves by the analy- 
sis of the scattering within a multiplet, characterized 
by the quantum number L. Initially this multiplet is 
(2L + 1)— fold degenerate. The impurities lower the rota- 
tional symmetry and lift the degeneracy. The transitions 
between the multiplets with different L are neglected. 
The matrix elements of the Hamiltonian are the projec- 
tion of the impurity 5— functions onto the multiplet, and 
we have 



N imt 

E 

=i 



Ui ro 2 Y£ m (n t )Y Ln (^) 



(48) 



In this section, in order to obtain a wider picture, 
we allow a variation of U{. Specifically, we numeri- 
cally consider three possibilities : i) the former case, 
all Ui = Mo, ii) Ui — ±uo and = 0, and iii) Ui is 
a normally distributed (ND) variable with variation u , 
cxp(ikuj) = cxp(— k 2 u 2 /2). 

Simple estimates show that the average off-diagonal 
matrix clement of Ti. scales in the limit of large Ni mp as 



(49) 



In this limit one may also regard different TL m ,n to be 
almost statistically independent quantities. Due to the 
known property of large random matrices, it follows then 
that the width of the impurity-induced band will be of 
order of upr^ 2 ^/Nj mp L. This estimate corresponds to 
the above one, eq. (p9|), and is validated by the numeri- 
cal calculations. At the same time, some features of our 
problem prevent us from an identification of Ti with one 
of the classes of random matrices. Hence the usual ex- 
pectations from the random matrix theory might not be 
fully applicable to our case. We discuss these points in 
the next subsections. 



A. Lifshits's theorem 

There exists a particular property, satisfied when the 
degeneracy (2L+1) is larger than the number of impurity 



centers, Ni mp . In this case at least 2L + 1 



Nimp levels 



remain degenerate for any choice of Ui,Qi. This state- 
ment is a general one for the projecti ve p erturbations 
and dates back to the work by Lifshits |19|. It is known 
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particularly for the impurity scattering of the electrons 
in a magnetic field. 15 lj] For a sphere, this statement 
can be proved as follows. 

Let us write Yr, m (^i) = {A m ) an d introduce a scalar 
product for two points on the sphere as 

or i i 



We introduce also a real- valued matrix M. with elements 
Mij = Uiro 2 {i\j). (50) 

Then TrH = J2 im u^ 2 {m\i) {i\m) = TyM, and Trft 2 = 

E l imn^ u j r o 4 ( n l i )(*l TO >( TO b'>(j|"> = TrA^ 2 , and gener- 
ally 

TrH k = TrM fc . (51) 
After this observation it is straightforward to show that 



det{E -H) = exp [{2L 
= exp [(2L 
= E 2L+1 



l)lnE + Ti-ln(l-H/E)} 
l)ln£ + Trm(l- M/E)] 
n 'det{E-M). (52) 



In the last line the pole in E at 2L + 1 < Ni mp is com- 
pensated by the necessary degeneracy of the matrix M. 
in this case. 

Therefore a singularity in the S-matrix of perturbation 
theory is identified with the incomplete raising of degen- 
eracy. 

Another specific feature of our problem is the follow- 
ing. We see that non-zero eigenvalues of ( f48| ) coincide 
with those of (|5^) which is a real symmetric matrix (for 
Ui = Uq) with a randomness in its elements. Then one 
should seemingly classify M to Gaussian orthogonal en- 
semble (GOE). However, the elements of M. are not sta- 
tistically independent. Indeed, it can be easily shown, 
that averaging over positions of impurities gives 



M 2 



(2L- 



l)u 2 



(4tt) : 



Q 2 



while, e.g., the triple combination with the cyclic se- 
quence of indices yields 



MijMjkMki 



V2L + 1 



7^0. 



This fact of nonvanishing long-range correlations makes 
our problem distinct from the usually considered ones. 
Being small, these correlations however should cause only 
minor deviations from the GOE predictions. 

For the case, when the amplitudes Ui are sign-reversal, 
M. is not symmetric and one expects larger deviations 
from the known picture. 



VimgUQ. As a result, the DOS 



B. numerical results 

We performed the numerical diagonalization of the ma- 
trices of the form (|48|), with the positions fij randomly 
distributed on the sphere and the amplitudes itj chosen 
according to one of the above ensembles. The calcula- 
tions were done with the use of standard EISPACK rou- 
tines for 10,000 random realizations of fli,Ui for each 

E, ^imp- 

For given L, Ni mp , we determine the density of states 
p(E) (DOS), averaging over the realizations. The density 
is normalized as p{E) dE — 2L + 1 . Some typical 
curves for DOS are shown in the Fig. |[ We observe 
following qualitative features : 

i) For identical impurities, Ui = uq, we have a smooth 
asymmetric DOS at A > 1, when the degeneracy is re- 
moved. At the point X = 1, DOS exhibits a mild sin- 
gularity, roughly of the law In 2 E. When X < 1, a part 
of the multiplet is unsplit, and the split-off states are 
separated from the position of this 6— function. For a 
convenience of presentation of Fig. |^a, we shifted the en- 
ergies by TtH/(2L + 1) 

is centered for each X, i.e. /_ E p{E) dE = 

ii) For the "dichotomic" impurities, = ±uq, ul = 0, 
the DOS is symmetric, showing the same features as in 
the case (i). 

iii) For the case of the normally distributed m, the prop- 
erties of DOS arc similar to the case (ii), except for the 
region X < 1. In this latter region, p(E) shows a mild 
singularity at E = 0, accompanied by S— function contri- 
bution of unsplit states. 

In order to test the accuracy of calculations, we fixed 
the window for the expected Nq = 2L + 1 — Ni„ lp unsplit 
states to the size 10~ 6 of estimated band-width. The 
average number of energy levels, found in this window 
was No + 5N with 8N < 1. In the worst case of ND 
amplitudes m and X < 1 we had SN < 10~ 3 . The small 
values of SN allowed us to separate the smooth contri- 
bution from the S— core. 

As a next step, we unfolded the spectrum, {Ei} — > {e;} 
by defining = J_ * p(E) dE. A certain subtlety in this 
procedure should be described here. If the degeneracy 
of the multiplet is completely removed, then the aver- 
age spacing Si = — £; between the nearest neigh- 
boring £j is s? = 1. We examine also the statistics of 
the spectrum for the partially degenerate level when the 
DOS contains a smooth part and a singular contribution, 
p(E) = p sm ooth(E) + N S(E). To describe this case, we 
exclude the states in the S— core from our consideration 
and write for the split levels £j = J_* p S mooth{E) dE. 

Obviously, upon doing this, we have again (ej+i — £j) = 
1 , now with i = l,..., Ni mp — 1 • 

We studied the energy level spacing distribution func- 
tion P(s) with the use of the unfolded spectrum. We 
remind, that usually in analyzing the spectrum of ran- 
dom ensembles, one finds either the Wigner-Dyson (WD) 
law, P(s) ~ sexp(— s 2 ), established for the GOE en- 
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semble, or the Poissonian law, P(s) ~ e~ s . Besides, 
it was shown that in the vicinity of the metal-insulator 
Anderson transition (MIT) an intermediate distribution 
may take place, which interpolates between the WD and 
Poissonian dependence. [^0) Particularly, one observes a 
"semi-Poisson" distribution, P(s) = 4se~ 2s in three di- 
mensions (D) at MIT. J2l| Similar intermediate P(s) oc- 
cur also in the studies of the statistics of the lowest Lan- 
dau level in the presence of disorder |22| and in some 
other models. Q 

It turned out, that our data are well described (with 
some predictable exclusions, discussed below) by the fol- 
lowing formula 

P(s) = cs exp(-as fc ). (53) 

Since the conditions J P(s)ds = J s P(s)ds = 1 uniquely 
determine the parameters a, c in the form 

a =(r[3/6]/r[2/6])\ (54) 
c = bT[3/b} 2 T[2/b}- 3 , 

the eq. (|5^) depends, in fact, only on b. This parame- 
ter was determined by fitting the numerically obtained 
curves both for P(s) and In P(s). Remarkably, in those 
cases, when a good visual agreement of the data with a 
fitting curve was achieved, the values of 6, determined 
both ways, essentially coincided. Some of the results are 
shown in Fig. |, |, |, [| 

The values of b are summarized in the Table Q. In 
cases of good visual agreement, the corresponding error 
in b was ±0.01, estimated from the different fitting pro- 
cedures used. From this table, we see that the values of b 
in the investigated range of L, Ni mp depend only on the 
ratio of these quantities. The obtained 6s are essentially 
the same in each of two domains X > 1 and X < 1, 
with somewhat lower 6s at X < 1. The worse agreement 
with the fit was obtained in the cases, when the regular 
part of DOS, p sm ooth(E) showed a mild singularity, Fig. 
||. In these cases we checked the procedure of unfolding 
the spectrum, which included the cubic spline interpola- 
tion. The range of E was divided, so that the position of 
singularity corresponded to a free boundary of a spline. 
The resulting spline curve fitted p sm ooth well, showing 
no oscillations , however, it did not improve the qual- 
ity of the fit of P(s) by eq. (|53|). Particularly, the tail 
of P(s) in these cases showed roughly exponential decay 
In P(s) ~ —s at larger s and a small non- vanishing value 
of P(s = 0), which can be expected for the systems with 
a localization of the states. The results of the fit of this 
tail by the exponential law for the normally distributed 
Ui are shown in Fig. ||. It is seen here, that P(s) at large 
s decreases faster with the increase of disorder X. 

The linear tail in In P(s) at large s for the singular DOS 
can be, in principle, expected. One can argue that in this 
case the unfolding procedure itself is not particularly use- 
ful. Indeed, this procedure assumes that different parts 
of the spectrum can be treated on equal grounds. This 



assertion can be tolerated for unsingular DOS, while the 
presence of a singularity in p sm ooth explicitly breaks the 
premise for the unfolding, Ei — ► e,. As a result, one 
anticipates that the part of the spectrum at the center 
of the band shows the intermediate statistics, which is 
different from the Poissonian one, stemming from the lo- 
calized states at the tails of the band (see, e.g., Ref. plfl). 

Note, that the considered cases of dichotomic and ND 
amplitudes lead to the same value of 6 = 1.76 ± 0.02 at 
X > 3. By analogy with the integer quantum Hall ef- 
fect (IQHE) |24j one can argue that the limit of large X 
corresponds to the white noise distribution of impurity 
potential ( X is associated with the number of impurities 
per one flux quantum). Our findings suggest that 6 =/= 2 
in this limit. 

Next, we performed the analysis of the number vari- 
ance, 

^(e)=W^)-(N{e)) 2 , (55) 

which measures the fluctuation of the number of levels 
N(e) in a band of unfolded spectrum of width e. The 
number variance is generally believed to be more sensi- 
tive than P(s) to the change of ensemble statistics. One 
has X 2 (e) = £ for the Poisson sequence, while GOE pre- 
diction is X 2 (e) ~ 27r~ 2 lne at large e. For the inter- 
mediate statistics at MIT, earlier theoretical arguments 
p5| , ^6| suggested that S 2 (e) oc e 2 ~ b , later it was agreed 
that X 2 (e) - xe with X < 1- 

Since we deal with finite size matrices, the range of e 
is restricted from above by (2L + 1). Actually the finite- 
size effects are felt at e ~ 0.4(2L + 1), where one finds a 
maximum of X 2 . The behavior of the number variance 
at lower e is shown in Fig. |[ 

From this figure one sees that in the case (i) of identi- 
cal impurities X 2 follows the GOE prediction except for 
the value X = 1. In cases (ii) and (iii), when are al- 
lowed to vary, the dependence X 2 (e) goes visibly above 
the GOE law. One notices that in the case (iii) of ND ui 
the obtained X 2 is considerably higher than in the case 
(ii) of dichotomic amplitudes. The highest curve for the 
number variance is obtained again at X = 1, both for 
the case (ii) and (iii). 

The pronounced variation of X 2 with X in two latter 
cases should be contrasted with our results for P(s), Ta- 
ble |. Indeed, we see that the the curves for X 2 may 
almost coincide for different 6's (at X < 1 and X > 1), 
while they may be notably different for almost identical 
6's. Particularly, one finds almost linear dependence at 
e > 3 , X 2 (e) ~ x£ + const with \ < !• The slopes x var y 
in the range 0.06 -j- 0.14 and 0.1 -j- 0.4 for the dichotomic 
and ND amplitudes, respectively. 

Making this observation we try to compare our system 
with the problems of metal-insulator Anderson transition 
and of impurity-split lowest Landau level. We mentioned 
the earlier theoretical expectation X 2 (e) oc e 2 ~ b above. 
This possibility, numerically validated for the 3D MIT 
pPf is apparently not favored by our data. Other studies 



7 



support the viewpoint that S 2 (e) oc xN with the level 
compressibility x < 1 ( see Pl| for a review). A rela- 
tion 2\ = 1 — Di/d connects J29| , [32| the quantity x with 
the multifractal exponent D2, characterizing the spatial 
extent of the wave functions in the flat geometry of d 
dimensions, J dr\^\ 4 oc V~° 2 / d with V a volume of the 
system. For the IQHE various authors pro vided the val- 
ues D 2 = 1.4 -7- 1.6 and x = 0.1 -7- 0.15. |33^|5 

The values of x extracted from our data are roughly 
consistent with these results. A pronounced dependence 
of x on X> however, requires further understanding. We 
suggest here that 

i) either our results obtained for e < 15 do not describe 
the true asymptote of £ 2 , (cf. @|l[) 

ii) or the multifractal exponent D2 varies with X. 

It would be interesting to check the last possibility by 
investigating the spatial character of the wave functions 
in our model. This work is however beyond the scope of 
the present study. 



IV. DISCUSSION AND CONCLUSIONS 

We considered the electron gas on the sphere, moving 
in the potential field of point-like scatterers. Without a 
disorder, one observes a series of effects, associated with 
the coherent motion of the electrons in the spherical ge- 
ometry. The high symmetry of a problem leads to the 
large degeneracy of energy spectrum in the ideal case. 
The disorder removes the rotational symmetry, prompt- 
ing one to expect the disappearance of the coherence ef- 
fects. 

Our consideration shows that there exist a definite 
range of disorder parameters, where the sphere remains 
sphere at least partly. Particularly we demonstrate that 
the phenomenon of incomplete raising of degeneracy may 
take place for sufficiently small radii. In the case of rel- 
atively small disorder the main effect is the splitting of 
the rotational multiplet, while the transitions between 
the different multiplets can be ignored. Numerically in- 
vestigating this latter case, we studied the energy level 
statistics in the impurity-induced band. We found that 
the energy spacing distribution function follows a modi- 
fied Wigner-Dyson law, P(s) ~ sexp(— s b ) with b ~ 1.76 
in the white noise limit. Our data suggest that the pa- 
rameter b does not determine the behavior of the energy 
levels' number variance E 2 . The latter quantity reveals 
various possibilities for almost the same values of b, de- 
pending strongly on the chosen model of disorder. 
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TABLE I. The values of the best fit to P(s) by the for- 
mula P(s) = cs exp(— as b ). Parameters a,c depend on b as 
described in text. 



(2L + 1) N mp X b( Ul = 1) b(uj = ±1) b(uj = ND) 



61 


300 


4.92 


1.88 


1.74 


1.78 


41 


200 


4.88 


1.89 


1.76 


1.76 


41 


120 


2.93 


1.89 


1.75 


1.75 


41 


50 


1.22 


1.86 


1.73 


1.62 a 


61 


61 


1 


1.64 1 


1.68 


1.34 a 


61 


50 


0.82 


1.80 


1.68 


1.45" 


81 


50 


0.62 


1.82 


1.68 


1.43 a 



a worse visual agreement with P(s), linear tail of In P(s) at 
s>2 



TABLE II. The values of the best fit to In P(s) by the for- 
mula In P(s) = — oi — a2(s — 2) at 2 < s < 5. ND amplitudes. 



(2L + 1) 


N imp 


X 


ai 


a-2 


41 


50 


1.22 


1.75(6) 


2.75(4) 


61 


61 


1 


1.81(2) 


2.20(1) 


61 


50 


0.82 


1.97(5) 


2.12(3) 


81 


50 


0.62 


2.15(8) 


1.82(3) 
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FIG. 3. The average density of states as a function of en- 
ergy, measured in units discussed in the text. The data for 
identical (a), sign-reversal (b) and ND amplitudes (c) arc 
shown. The values of X are X = 4.92 (□), 1 (o), and 0.62 (A), 
For convenience of presentation of data with different L, the 
shown DOS is scaled as J p{E) dE — 1 here. The lines are 
guides to the eye, every second data point is shown. The verti- 
cal line on the plot (a) denotes the position of the 8— function 
for X = 0.62. 
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FIG. 4. Energy spacings distribution for L = 61, 
: 300, identical impurities; b = 1.88 . 




FIG. 5. Energy spacings distribution for L 
Nim P = 50, identical impurities; b = 1.82 . 
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FIG. 7. 



Energy spacings distribution for L 



20, 



FIG. 6. 



Energy spacings distribution for L = 30, 



50, "dichotomic" impurities; b— 1.68 
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200, ND impurities; b = 1.76 
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FIG. 8. Number variance for the amplitudes with (a) 
fixed value, (b) "dichotomic" and (c) normal distributions. 
The values of X are : 4.92(D), 4.88 (o), 2.93(A), 1.22 (v), 
1.0 (O), 0.82 (+), 0.62 (x). The GOE prediction for the num- 
ber variance is shown by a solid line. 



11 



